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I .  INTRODUCTION 


fHECEfllflG  PAGE  BUNK  -  NOT  J3XMBD 


Compressive  heating  caused  by  the  rapid  collapse  of  an  air  layer 
adjacent  to  a  high  explosive  material  has  been  implicated  as  a  mechanism 
of  ignition  under  artillery  setback  conditions.  This  conclusion  has 
been  reached  as  a  result  of  activator  experiments  conducted  at  the  Bal¬ 
listic  Research  Laboratory  (BRL)1.  The  activator  is  an  experimental  appa¬ 
ratus  in  which  an  air  gap  adjacent  to  a  high  explosive  sample  is  rapidly 
compressed  by  a  piston.  In  support  of  these  experiments,  two  one-dimensional 
models  have  been  developed  to  trace  the  thermal  history  of  the  explosive 
sample,  the  air  gap  and  the  piston. 

The  apparatus,  illustrated  schematically  in  Figure  1,  consists  of  a 
heavy  confinement  cylinder  with  a  12.7mm  bore  enclosing  the  explosive 
sample, a  driving  piston  and  a  backup  piston.  The  driving  piston  is 


Figure  1.  Activator  Schematic 


1 Taylor ,  B. ,  C. ,  Starkenberg ,  J. ,  and  Ervin ,  L.  H. ;  "An  Experimental 
Investigation  of  Composition  B  Ignition  under  Artillery  Setback 
Conditions  3  BRL  Report  to  be  published. 
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set  in  motion  by  a  larger  piston  which  is  in  turn  driven  by  a  propellant 
burned  in  a  low  pressure  chamber  or  breech.  The  backup  piston  rests 
against  a  rigid  stop.  An  air  gap  is  set  between  the  explosive  and  the 
driving  piston  and  an  additional  free  run  gap  may  be  allowed  between  the 
driving  and  the  large  piston.  The  explosive  sample  is  heated  as  the 
air  gap  is  compressed.  The  sum  of  the  air  gap  thickness  and  the  free 
run  is  referred  to  as  the  total  run.  In  order  to  control  the  gap  closure 
rate,  shear  pins  are  used  to  allow  the  breech  to  pressurize  before  the 
large  piston  begins  to  accelerate.  Thus  the  pressurization  rate  and 
peak  pressure  depend  on  propellant  charge  design,  total  run  and  shear  pin 
strength.  It  is  further  possible  to  pressurize  or  evacuate  the  air  gap 
and  to  use  various  insulating  materials  at  the  face  of  the  driving  piston. 

The  experimental  evidence  to-date  indicates  that  the  sensitivity  of 
the  explosive  to  violent  reaction  is  governed  by  a  number  of  factors. 
Specifically  these  are: 

1.  total  run 

2.  peak  pressure 

3.  pressurization  rate 

4.  initial  air  gap  thickness 

5.  initial  air  pressure 

6.  piston  thermal  conductivity 

7.  air  leakage 

8.  state  of  the  explosive  surface 

In  the  experiments  conducted  thus  far  the  first  three  of  the  above  were 
not  varied  independently.  The  list  does  not  exhaust  all  possible 
influencing  factors.  The  present  study  addresses  all  but  the  last 
of  these. 

A  preliminary  and  final  model  have  been  constructed.  For  the  pre¬ 
liminary  model  the  air  gap  is  assumed  to  be  instantaneously  compressed 
to  its  peak  pressure  and  temperature.  It  is  then  possible  to  compute 
the  temperature  of  the  explosive-air  interface,  which  is  the  parameter 
used  for  ranking  the  sensitivity  of  various  configurations  in  this  study. 
This  model  corresponds  to  an  intuitive  approach  to  the  problem  referred 
to  as  adiabatic  compression  which  has  been  widely  used  to  describe  comp¬ 
ressive  heating.  This  concept  is  inappropriate  when  applied  to  the  entire 
problem  since  it  identifies  only  a  limiting  case  which  is  incomplete  and 
can  be  misleading.  In  order  to  more  completely  describe  the  situation 
a  model  in  which  the  air  gap  is  compressed  at  a  finite  rate  is  required. 

In  this  more  sophisticated  model,  the  heat  equation  is  solved  numerically 
for  the  temperature  distribution  in  the  three  layers  as  a  function  of  time. 
Compressive  heating  is  introduced  through  a  source  term  active  in  the  air 
layer.  The  peak  explosive-air  interface  temperature  is  extracted  from 
these  computations.  An  implicit  finite  difference  scheme  has  been 
employed  because  of  its  flexibility  and  applicability  to  parabolic  prob¬ 
lems.  The  adiabatic  compression  model  may  be  regarded  as  the  limiting 
case  for  infinite  pressurization  rates. 
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The  objective  of  the  modeling  is  to  treat  the  ignition  problem 
by  examining  the  effects  of  stimulus  variations  on  sensitivity  in  order 
to  verify  the  compressive  heating  hypothesis  and  to  qualitatively  predict 
activator  behavior.  A  qualitative  approach  is  necessary  because  the 
models  are  limited  by  their  one-dimensional  property  while  events  occurring 
in  the  activator  and  in  artillery  ammunition  are  multidimensional  in 
nature.  Steps  have  been  taken  to  account  for  some  of  the  multidimensional 
aspects  of  the  problem. 


II.  MODELING 

A.  Adiabatic  Compression  Model 

As  a  first  step  consider  the  situation  that  arises  when  the  air  gap 
is  instantaneously  compressed  from  a  known  initial  pressure  to  a  known 
final  pressure.  Assuming  that  the  air  is  an  ideal  polytropic  gas  results 
in  the  following  expression  for  final  air  temperature,  density  and  gap 
thickness . 


Y-l 


The  temperature  profile  immediately  after  compression  is  in  the  form 
of  a  step  with  ambient_temperature  persisting  in  the  explosive  and  the 
elevated  temperature,  Tp  in  the  air.  In  this  case  the  highest  temperature 

arising  in  the  explosive  is  the  explosive-air  interface  temperature.  This 
can  be  calculated  using  a  solution  to  the  heat  equation  for  two  semi¬ 
infinite  layers  instantaneously  brought  into  contact  following  Carslaw 
and  Jaeger2.  The  resulting  expression  for  interface  temperature  is 
independent  of  time. 


Y-l 


2  Carslaw,  H.  S.  and  Jaeger,  J.  C.;  Conduction  of  Heat  ■ in  Solids, 
2nd  ed. ,  Oxford  University  Press,  1959,  pp  87  -  88. 
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where 


The  following  expression  for  the  ratio  of  the  increase  in  interface 
temperature  to  the  increase  in  air  temperature  may  be  generated. 


This  is  referred  to  as  the  heating  efficiency  of  the  gas  with  respect  to 
the  explosive.  Furthermore,  a  is  a  function  of  the  final  temperature  and 
density  of  the  air  and  thus  of  the  compression  ratio,  p^/p^,  and  the  initial 

state  of  the  air  before  compression.  It  is  possible,  then,  to  predict  the 
interface  temperature  knowing  the  compression  ratio  and  the  properties 
of  the  explosive  and  the  air. 


This  model  excludes  effects  of  pressurization  rate,  gap  thickness, 
piston  insulating  capacity,  and  mass  addition  or  leakage;  these  are 
studied  by  means  of  the  more  sophisticated  model.  A  slightly  more  detailed 
approach  which  utilizes  a  solution  to  the  heat  equation  for  all  three  layers 
could  be  used.  In  this  event  the  same  initial  interface  temperature 
should  be  predicted  since  for  sufficiently  early  times  after  compression 
even  a  very  thin  air  layer  will  act  as  though  infinite.  At  later  times 
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the  effect  of  the  third  layer  would  be  to  reduce  the  interface  temperature. 
Thus  the  same  maximum  interface  temperature  is  predicted  and  no  additional 
information  is  generated  through  this  approach. 


B.  Finite  Compression-Rate  Model 


1.  Governing  Equations.  The  heat 
layers  consisting  of  the  explosive,  the 
trated  in  Figure  2.  It  takes  a  slightly 
layers.  An  Eulerian  representation  has 

equation  is  applied  in  three 
air  gap  and  the  piston,  illus- 
different  form  in  each  of  the 
been  used. 
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is  the  particle  derivative  and  J  represents  heating  due  to  chemical  re¬ 
action  in  the  explosive.  The  equations  may  be  generalized  as  follows: 


DT 

-  c 

iU| 

(7  nv 

®1‘- 
Dt  J 

Dt 

c  1 

P  |_3x 

'  3x  ' 

P 

where 

{  (0,1) 

explosive 

CV 

a2}  = 

{  (1,0) 

air 

(  (0,0) 

piston 

It  is  assumed  that  dissipation  may  be  neglected  in  the  air  layer.  The 
following  additional  assumptions  are  used  in  order  to  compute  conditions 
in  the  air  layer: 

1.  air  pressure  is  uniform  throughout  the  gap  and  varies  with 
time  in  a  prescribed  manner 

2.  the  ideal  gas  equation  of  state  applies 

3.  velocity  varies  linearly  in  the  air  layer 

4.  thermal  conductivity  is  constant  in  the  explosive  and  piston 
and  is  a  given  function  of  temperature  in  the  air  layer 


The  first  assumption  implies  that  the  pressurization  is  sufficiently 
slow  to  allow  wave  propagation  through  the  entire  air  layer  during  a  small 
increase  in  pressure.  The  consequences  of  assuming  a  linear  velocity/ 
profile  are  discussed  in  the  appendix.  Note  that  constant  density 
and  specific  heat  are  not  implied.  With  these  assumptions  the  general 
equation  becomes 


In  order  to  handle  the  compression  of  the  air  layer  an  expression 
for  the  total  mass  per  unit  area  of  the  air  in  the  gap  at  any  time  is  used 


To  account  for  leakage  or  mass  addition  during  piston  closure,  m  may  be 
regarded  as  a  function  of  time. 


The  ideal  gas  equation  of  state  is  used  to  relate  density  to  pressure 
and  temperature. 


Since  pressure  is  assumed  uniform 


This  expression  relates  the  pressure,  temperature  profile  and  gap  thick¬ 
ness  at  any  instant. 

The  particle  velocity,  u,  vanishes  in  the  explosive,  varies  linearly 
in  the  air  and  is  uniform  in  the  piston  at  any  time. 

explosive 
air 

piston 


The  thermal  conductivity  of  air  is  given  by 


K  . 


air 


(3) 


The  values  of  and  C 2  were  obtained  by  fitting  to  data  found  in 

reference  3  and  yield  an  accuracy  of  5  percent  between  300K  and  400K  and 
2  percent  between  400K  and  1500K.  This  covers  the  entire  range  of  gas 
temperatures  encountered.  The  effect  of  pressure  has  been  neglected. 

2.  Boundary  Conditions:  At  interior  boundaries  (those  between  layers) 
temperature  and  heat  flux  are  continuous. 

3  National  Bureau  of  Standards,  Tables  of  Thermal  Properties  of  Gases, 
Circular  546,  1955. 
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where  the  superscript  identifies  the  layer.  At  each  exterior  boundary 
(extreme  left  or  right)  either  temperature  or  heat  flux  may  be  specified 
as  a  function  of  time. 


left: 
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right: 
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3. 

Boundary  Motion.  The 

explosive-air 

interface  is  assumed 

to  lie 

at  a  fixed  point  so  that 


0  =  o/2> 


o 


T  i 

The  air-piston  interface  moves  at  the  piston  velocity 

0  (2)  =  o =  an 

r  Jt  p 

and  the  left  and  right  exterior  boundaries  advance  through  the  explosive 
and  the  piston  respectively  in  order  to  move  with  the  heated  region  in 
each  layer.  This  spreads  the  temperature  variation  over  the  portion  of 
the  layer  considered  at  any  time  and  leads  to  more  accurate  computations. 
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The  boundary  positions  are  then  given  by 


-M)  .  -(2) 

r  "  xl 


i(2)  =  ^ 

r  x* 


(3)  ,  5  *  J _ E_ 


,1/2 


x'-'  *  6  +  5| 

r 


lP  C 
1  p  pr 


(4) 


This  procedure  insures  that  the  explosive  and  piston  computational 
layers  are  maintained  at  the  appropriate  thickness  for  accurate  numerical 
computation. 


4.  Nondimensional  Forms.  Consider  the  following  coordinate 
and  variable  transformations: 
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When  introduced  into  equation  2  there  results  a  nondimensional  form  of 
the  heat  equation. 


3T 

3t 


p  cp  F 


3  ,  3T  ,  .  F  U  3T  al  '  2  T 

(<_  )  +  -  +  —  p  ♦  —  J 

F  3X  pCp  Cp 


2  3x  1  3x 


where 


u  dt 
o 


•  dF  Ur-U«  „  ,, 

F“  3t - : —  ■  ur  -  °t 


The  pressure  relation  can  also  be  written  in  nondimensional  form 

dx 


f 

JO 


m  =  p  F  J  ip-  =  pFV 


where 


R  .  T 

_  air  o  - 

m  =  — - -  m 

P  L 
ro  o 


The  expression  is  solved  for  F  with  pressure  and  mass  given, 


and  for  F 


(5) 


where 


(  m  _  1_ 
^  m  y 


W 

m  cp 


32T  3£  3T  .d^ 

„  2  +  3x  3x  T 

3x 


(6) 


Mass  variation  during  gap  closure  enters  the  model  only  through  equations 
5  and  6.  This  implies  the  assumption  that  mass  is  varied  throughout  the 
air  layer  at  the  local  temperature  such  that  only  the  volume  is  affected. 
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The  nondimensional  form  for  air  thermal  conductivity  is 


where 
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air  1  2 


T  1/2  C. 
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The  interior  boundary  conditions  are  simply 
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right: 

The  nondimensional  forms  of  the  boundary  velocities  are 
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5.  Finite  Difference  Equations.  The  computational  grid  is  shown  in 
Figure  3”  The  coordinates  of  the  grid  points  are  given  by 
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max 


1 <i <in 
- max 
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l« - REGION  n  - ► 


Figure  3.  Computational  Grid 


The  grid  includes  fictitious  points  on  either  side  of  each  region  for 
imposing  conditions  on  derivatives,  such  that  the  boundary  points  of  each 
region  correspond  to  i=2  and  i=imax-l.  The  pattern  for  the  computation  of 

derivatives  is  illustrated  in  Figure  4.  This  leads  to  an  implicit  finite 
difference  scheme  for  the  solution  of  the  equation.  The  finite  difference 
form  of  the  equation  is  applied  at  interior  points. 
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Figure  4.  Derivative  Pattern 


The  superscript  £  refers  to  values  at  the  current  time  step  and  the 
superscript  k  to  values  at  the  previous  time  step.  The  interior  bound¬ 
ary  conditions  are  written 
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The  total  number  of  equations  and  boundary  conditions  is  equal  to  the 
total  number  of  unknowns. 
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max 
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max 


Solutions  at  each  time  step  are  generated  by  computing  F  and  F  from 
equations  5  and  6  and  then  solving  the  set  of  simultaneous  equations 
represented  by  equations  7,  8,  and  9.  The  finite  difference  equations 
including  boundary  conditions  must  be  solved  in  all  three  regions 
simultaneously.  The  solution  requires  the  inversion  of  a  band  matrix  with 
3  upper  and  2  lower  codiagonals.  This  is  accomplished  by  use  of  a  suitable 
matrix  inversion  routine.  The  procedure  is  repeated  in  iterative  fash¬ 
ion  until  an  arbitrarily  specified  degree  of  accuracy  is  attained. 

III.  RESULTS  AND  DISCUSSION 


A.  General 


The  material  constants  for  the  explosive,  the  air  and  the  piston  are 
summarized  in  Table  I.  For  the  explosive,  values  for  Composition-B  (Comp-B) 
are  given  and  for  the  piston  those  for  polyethylene  and  steel  are  given. 

In  addition  to  values  for  air,  those  for  argon,  carbon  dioxide  and  hydrogen 
are  included.  The  values  of  specific  heat,  obtained  from  reference  3, 
correspond  to  atmospheric  pressure.  The  value  for  air  applies  at  1300K 
while  those  for  argon,  carbon  dioxide,  and  hydrogen  apply  at  300K.  The 
initial  temperature  of  the  system  is  always  taken  to  be  300K. 


Table  I 

.  Material 

Constants 

pi 

c 

P 

K 

C1 

C2 

kg/m3 

J/kgK 

W/mK 

CM 

\ 

to 

US 

e 

W/mK 

Comp-B 

1.69xl03 

1.393xl03 

.2623 

0. 

.2623 

Air 

1.18 

1.197xl03 

_ 

2.901xl0"3 

2.536x10"' 

Argon 

1.62 

5. 23xl02 

- 

1.743xl0'3 

1.256x10"' 

Carbon  Dioxide 

1.80 

8. 23xl02 

_ 

3.236xl0'3 

3.616x10"' 

Hydrogen 

8. 15xl0-2 

1 ,42xl04 

- 

1.928xl0'2 

1.523xl0'] 

Polyethylene 

9. 30xl02 

2.301xl03 

.4184 

0. 

-.4184 

Steel 

7.87xl03 

4.372xl02 

79.5 

0. 

-  79.5 
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When  the  time  scale  does  not  vary,  the  magnitude  of  the  peak  temperature 
attained  at  the  explosive-air  interface  is  indicative  of  the  likelihood 
that  a  given  configuration  will  lead  to  ignition.  In  the  subsequent 
discussion,  variations  which  increase  this  temperature  are  referred  to  as 
sensitizing  while  those  that  reduce  it  are  referred  to  as  desensitizing. 

In  the  case  of  the  adiabatic  compression  model,  this  is  the  temperature 
achieved  immediately  upon  piston  closure.  The  final  to  initial  tem¬ 
perature  and  density  ratios  for  y  =  1.4  are  plotted  as  functions  of 
compression  ratio  in  Figure  5a  and  the  resultant  interface  temperatures 
are  plotted  in  Figure  5b.  Since  the  model  assumes  a  constant  thermal 
conductivity  for  air  and  this  actually  varies,  upper  and  lower  limits 
for  the  interface  temperature  based  on  the  maximum  and  minimum  values  of 
k  have  been  computed.  The  maximum  value  occurs  at  the  final  air  tem¬ 
perature  and  the  minimum  occurs  at  the  interface  temperature.  The  results 
indicate  that  increasing  peak  pressure  is  sensitizing  while  increasing 
the  initial  pressure  is  desensitizing.  The  latter  conclusion  is  in  conflict 
with  the  experiments.  This  indicates  that  the  processes  involved  are  not 
adiabatic.  Further,  this  model  does  not  predict  sufficiently  high  temp¬ 
eratures  to  explain  ignitions  observed  at  the  lower  compression  ratios 

It  is  estimated  that  temperatures  in  excess  of  650K  sustained  for  several 
hundred  microseconds  are  required  for  ignition  of  Comp-B  in  the  activator. 
With  initially  pressurized  air  gaps,  ignitions  were  observed  with  com¬ 
pression  ratios  as  low  as  600  for  which  the  adiabatic  compression  interface 
temperature  lies  between  425K  and  502K. 

For  the  finite  compression-rate  model  the  heat  flux  was  required 
to  vanish  at  the  exterior  boundaries.  The  motion  of  these  boundaries 
through  the  material  results  in  the  maintenance  of  a  temperature  very 
close  to  300K  at  the  left  and  right  extremes  throughout  the  computation. 
Twenty-three  computational  points  were  used  in  each  region.  The  peak  ^ 

temperature  was  found  to  be  quite  sensitive  to  the  time  step  and  At=2xl0 
was  required  for  accuracy.  A  constant  pressurization  rate  during  piston 
closure  is  assumed. 


t  <  t 


f 


t 
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Figure  5.  Adiabatic  Compression  Computations 
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A  similar  expression  is  used  when  there  is  mass  variation. 


m  = 


">i  IF  -l  )  -  +i] 

t„ 


m. 

1 


t  <  t. 


tf  <  t 


where  t^,  the  gap  closure  time,  is  300  ps.  This  pressure  rise  time  appears 
to  be  characteristic  of  the  activator  system  when  the  total  run  is  not 

4 

equal  to  zero.  Thus  some  1.5x10  steps  are  required  to  reach  peak  pressure. 
Typical  evolution  of  the  temperature  profile  during  gap  closure  is  shown 
in  Figure  6a  and  typical  interface  temperature  and  maximum  air  temperature 
histories  are  compared  in  Figure  6b.  The  computations  are  terminated  shortly 
after  the  interface  temperature  peaks.  It  can  be  noted  that  the  maximum 
air  temperature  is  significantly  higher  than  the  interface  temperature 
and  that  it  peaks  earlier. 

B.  Effects  of  Parameter  Variations. 

A  number  of  computations  were  made  to  assess  the  influence  of  variations 
which  were  shown  to  affect  sensitivity  in  the  activator  experiments. 
Specifically  peak  pressure  and  pressurization  rate,  initial  pressure, 
initial  gap  thickness,  piston  material,  air  leakage  during  gap  closure 
and  the  presence  of  alternate  gases  in  the  gap  were  considered.  The 
results  of  this  are  summarized  in  Tables  II  through  VII.  Case  B  is 
common  to  all  the  tables  and  the  fixed  conditions  for  each  table  correspond 
to  that  case. 

1 .  Peak  Pressure  and  Pressurization  Rate. 

When  a  fixed  value  of  the  gap  closure  time  is  selected  the  peak 
pressure  and  pressurization  rate  are  coupled.  The  question  arises  as  to 
which  of  these  most  strongly  influences  peak  temperature.  In  case  A, 
appearing  in  Table  II  peak  temperature  and  pressure  occur  simultaneously. 

In  cases  B  and  C  the  effect  of  heat  conduction  is  sufficient  to  cause 
the  temperature  to  peak  and  the  computation  to  be  terminated  before  the 
maximum  pressure  is  attained.  Thus,  for  pressurization  rates  in  excess 
of  340  GPa/s  and  for  other  conditions  corresponding  to  Table  II  peak 
temperature  is  a  function  of  pressurization  rate  exclusively.  The  effect 
of  increasing  pressurization  rate  is  sensitizing  as  the  experiments  have 
shown. 


Table  II.  Effects  of  Pressurization  Rate 


Case 

Peak 

Pressure 

Pressur¬ 

ization 

Rate 

Pressure 
at  Peak 
Temperature 

Peak 

Interface 

Temperature 

GPa 

GPa 

GPa 

K 

A 

.1 

3.4xl02 

.  10 

432 

B 

.3 

l.OxlO3 

.25 

512 

C 

1.0 

3.4xl03 

.47 

613 

2 .  Initial  Pressure. 

As  initial  pressure  is  increased  the  effect  is  at  first  sensitizing 
but  subsequently  desensitizing.  The  adiabatic  compression  model  predicted 
only  the  desensitization  while  in  the  experiments  only  the  sensitizing 
effect  was  observed.  The  transition  begins  approximately  where  peak 
pressure  is  achieved.  This  suggests  that  the  experiments  lie  in  the  regime 
in  which  sensitivity  is  determined  by  pressurization  rate. 


Table  III.  Effects  of  Initial  Pressure 


Initial 

Mass 
of  Air 

Pressure 
at  Peak 

Peak 

Interface 

Case 

Pressure 

in  Gap 

Temperature 

Temperature 

MPa 

kg/m2 

r?  a 

K 

B 

.10 

5 . 89xl0~3 

.25 

512 

D 

.15 

8 . 83xl0~3 

.30 

549 

E 

.20 

1. 18X10"2 

.30 

555 

F 

.30 

1.77xl0'2 

.30 

543 

3.  Inital  Gap  Thickness. 

Increasing  the  initial  gap  is  also  sensitizing,  as  observed 
experimentally.  Cases  D  and  H  may  be  compared  to  show  that  increasing 
the  initial  gap  thickness  is  still  sensitizing  when  the  mass  of  air  in 
the  gap  is  held  constant. 
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Table  IV.  Effects  of  Initial  Gap  Thickness 


Initial 

Mass  of 

Peak 

Gap 

Air  in 

Interface 

Case 

Thickness 

Gap 

Temperature 

m 

kg/m  2 

K 

G 

2.5xl0-3 

2. 94xl0-3 

428 

B 

5.0xl0-3 

5. 89xl0"3 

512 

H 

7 . 5xl0-3 

8.83xl0'3 

573 

4.  Piston  Material. 


Increasing  the  thermal  diffusivity  c  )  of  the  piston  by  sub¬ 

stituting  steel  for  polyethylene  is  only  mildly  desensitizing. 


Table  V.  Effects  of  Piston  Material 


Peak 

Thermal 

Interface 

Case 

Material 

Diffusivity 

Temperature 

Wm2/J 

K 

B 

Polyethylene 

1 . 96xl0-7 

512 

I 

Steel 

2.31xl0'5 

505 

5 .  Leakage. 

The  leakage  of  fifty  percent  of  the  air  during  gap  closure  is 
desensitizing.  The  experiments  indicate  that  using  a  tightly  fitting 
piston  increases  sensitivity. 


Table  VI.  Effects  of  Leakage 

Pressure  at  Peak 

Case  in^/rru  Peak  Temperature  Temperature 


B 

J 


1.0 

0.5 


GPa 

.19 

.25 


K 

512 

494 


6.  Alternate  Gases. 


Gases  other  than  air  have  been  substituted  by  some  experimenters4 
in  order  to  test  the  compressive  heating  hypothesis.  In  general  the 
results  have  been  inconclusive.  Here  computations  have  been  made  for  argon, 
carbon  dioxide  and  hydrogen.  Results  from  the  adiabatic  compression  model 
are  given  in  Figure  7,  where  only  the  lower  value  of  the  interface  temp¬ 
erature,  based  on  the  interface  value  of  k,  has  been  plotted.  This  model 
predicts  a  significant  difference  in  the  pressure  required  for  ignition 
with  different  gases.  The  results  of  finite  compression-rate  model 
computations  are  presented  in  order  of  increasing  sensitivity  in  Table  VII. 
The  order  of  sensitivity  is  the  same  for  both  models.  Argon  has  the 
highest  ratio  of  specific  heats  and  produces  the  hottest  gas  on  compression. 
However,  since  it  possesses  the  smallest  heating  efficiency  of  the  four 
gases  it  transfers  the  smallest  percentage  of  its  temperature  increase  to 
the  explosive  and  only  ranks  second  in  sensitivity.  The  high  heating 
efficiency  of  hydrogen  causes  the  adiabatic  compression  model  to  predict  a 
sensitivity  significantly  higher  than  for  the  other  gases  while  the  finite 
compression-rate  model  shows  only  a  very  slightly  higher  sensitivity.  Note 
that  the  high  thermal  diffusivity  of  hydrogen  causes  the  temperature  to 
peak  at  a  much  earlier  point  in  the  compression  when  the  gas  temperature 
is  lower. 


Table  VII.  Effects  of  Alternate  Gases 


Thermal 

Heating 

Pressure  at 

Peak 

Gas 

y 

Diffusivity 

Efficiency  Peak  Temperature 

Temperature 

Wm'YJ 

GPa 

K 

CO- 

1.31 

1.34xl0"5 

4. 77xl0*5 

.26 

490 

c 

Air 

1.40 

1.76xl0'5 

5. 70x10”“’ 

.25 

512 

Argon 

1.67 

2.08xl0"5 

2.42x10“’ 

.28 

628 

H2 

1.41 

1.57xl0“4 

3. 40xl0-4 

.09 

632 

Validity  of 

the  One-Dimensional  Models 

The 

time  to 

explosion  for 

the  activator 

experiments  is 

judged  to  be 

within  a  few  hundred  microseconds  of  the  time  at  which  peak  temperature 
is  reached.  By  using  kinetic  parameters  for  RI)X  it  is  possible  to  roughly 
estimate  that  sustained  temperatures  in  excess  of  650K  are  required  for 


4  Arthur  D.  Little ,  Inc.,  "Cavity  Standards  for  Cast-Loaded  Artillery 
Projectiles" ,  Revised  Final  Report  for  Pi-catinnj  Arsenal,  30  Mar  57. 
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ignition  on  this  time  scale.  The  adiabatic  compression  model  predicts 
the  highest  possible  values  of  the  explosive-air  interface  temperature. 
Reference  to  Figure  5a  indicates  that  sufficiently  high  temperatures  can 
only  be  produced  at  compression  ratios  higher  than  many  at  which  ignition 
was  observed.  In  addition,  at  finite  pressurization  rates  even  lower 
temperatures  are  predicted  and  in  no  case  can  observed  ignitions  be 
accounted  for.  The  situation  is  further  aggravated  by  the  fact  that 
air  leakage  in  the  activator  renders  the  environment  even  less  hostile. 

This  discrepancy  between  the  predictions  and  the  experiments  arises 
because  the  one-dimensional  planar  gap  model  is  not  strictly  valid.  Among 
the  multidimensional  effects  that  may  come  into  play  are: 

1 .  enhanced  energy  transport 

2.  rapid  pressurization 

3.  convergent  air  flow 

4.  dieseling 

In  order  to  critically  examine  how  these  influence  sensitivity  it  is  neces¬ 
sary  to  introduce  them  into  the  one-dimensional  scheme. 

1.  Enhanced  Energy  Transport.  One  of  the  assumptions  of  the  analysis 
is  that  thermal  conductivity  is  the  only  mechanism  of  energy  transport 
through  the  air  layer.  Tf  turbulence  develops  in  the  air  layer,  this 
will  enhance  energy  transport  and  may  serve  to  increase  the  interface 
temperature.  This  effect  may  be  readily  introduced  into  the  present 
models  by  assuming  that  the  air  possesses  an  effective  thermal 
conductivity  which  is  greater  than  the  actual  value.  This  is  reflected 
through  the  value  of  a.  in  equation  1  and  is  introduced  into  the  finite 
compression  rate  model  by  specifying  a  thermal  conductivity  multiplier, 
o,  defined  as  the  ratio  of  the  effective  to  the  actual  conductivity. 

Of  the  multidimensional  effects  considered,  only  enhanced  energy 
transport  serves  to  increase  the  theoretical  maximum  temperature  of  the 
adiabatic  compression  model  defined  by  equation  1.  Neither  pressurization 
rate  nor  degree  of  convergence  can  be  varied  so  as  to  produce  higher 
temperatures  than  this  during  finite  rate  compressions.  Included  in 
Figure  8  is  a  plot  of  theoretical  maximum  interface  temperature  versus 

3 

the  thermal  conductivity  multiplier  for  a  pressure  ratio  of  3x10  .  The 
temperature  varies  between  the  limits  of  initial  temperature,  T. ,  for  o.=0 

v-i  1 

and  final  air  temperature,  T^(p^./p^)  Y  ,  for  Thus  the  interface 

temperature  becomes  increasingly  sensitive  to  compression  ratio  with 
increasing  thermal  conductivity. 

When  the  finite  compression-rate  model  is  used  different  results  are 
obtained  depending  on  the  variation  of  energy  transport  enhancement  through 
the  air  layer.  The  expression  used  for  o  is 
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0  =  (W15  (1-X)S  +  1 

When  s=0  the  thermal  conductivity  is  augmented  by  the  factor  °max  uniformly 

throughout  the  air  layer.  For  larger  values  of  s  the  effect  is  con¬ 
centrated  nearer  to  the  explosive  surface  as  in  a  boundary  layer.  This  is 
illustrated  in  the  inset  in  Figure  8  where  the  explosive  lies  at  x=0  in 
the  air  layer.  If  there  are  no  boundary  effects  only  a  very  limited  in¬ 
crease  in  temperature  results  from  enhanced  energy  transport  as  illustrated 
in  Figure  8  for  s=0.  The  temperature  remained  below  650K  for  values  of 

0max  brought  106.  The  reason  that  this  occurs  is  more  clearly  illustrated 

in  Table  VIII.  As  a  is  increased,  peak  temperature  occurs  at  increasingly 

early  times  during  the  compression  since  the  thermal  conductivity  increases 
while  the  rate  of  heat  generation  from  pressurization  remains  unchanged. 

The  interface  temperature  is  then  limited  by  total  energy  availible  in 
pressurizing  the  air  to  increasingly  lower  pressures.  Thus,  in  the  absence 
of  boundary  effects,  enhanced  energy  transport  is  ineffective  in  increasing 
the  temperature  prediction  unless  the  conditions  of  adiabatic  compression 
are  simultaneously  approached.  This  may  be  accomplished  through  rapid  pres¬ 
surization  or  convergence.  With  concentrated  energy  transport  enhancement  (s>0) 
it  is  possible  to  predict  temperatures  high  enough  to  explain  ignition. 

with  values  of  a  less  than  ten  as  illustrated  in  Figure  8  for  s=l,2  and  5. 

max  6  ’ 


Table  VIII. 

Effects  of  Energy  Transport 

(s=0) 

a 

max 

Pressure  at 

Peak  Temperature 

Peak 

Temperature 

GPa 

K 

1 

.25 

512 

102 

.08 

571 

10 

.04 

607 

103 

.04 

609 

io4 

.04 

611 

105 

.04 

614 

io6 

.04 

615 

2.  Rapid  Pressurization. 

Since  uniformly  increased  thermal  conductivity  favorably  affects 
the  results  of  the  adiabatic  compression  model  but  not  the  results  of  its 
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nonadiabatic  counterpart,  it  is  of  interest  to  consider  an  effect  which 
causes  the  finite  compression-rate  model  to  approach  the  adiabatic 
compression  limit.  This  occurs  as  the  pressurization  rate  becomes 
arbitrarily  large  and  can  be  introduced  into  the  finite  compression-rate 
model  by  specifying  the  pressurization  rate  independent  of  the  peak 
pressure  when  making  the  computation.  This  has  the  desired  effect  on 
the  interface  temperature  prediction  as  illustrated  in  Figure  9  for 
a  =  omax  =  10.  These  computations  were  made  with  the  time  step  reduced 

in  inverse  proportion  to  the  pressurization  rate.  The  temperature  approaches 
the  adiabatic  limit  with  increasing  pressurization  rate.  The  computations 
could  not  be  carried  to  higher  pressurization  rates  without  a  significant 
loss  of  accuracy  since  an  increasingly  small  region  of  the  air  layer  is 
affected  by  heat  conduction.  Temperatures  sufficiently  high  to  explain 
ignition  are  obtained  by  increasing  the  pressurization  rate  by  a  factor 
less  than  ten. 


Figure  8.  Effect  of  Energy  Transport  on  Interface  Temperature 
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Figure  9.  Effect  of  Pressurization  Rate  on  Interface  Temperature 
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3.  Convergent  Air  Flow.  Reference  to  Figure  5a  indicates  that  for 
pressures  in  the  range  of  interest  a  gap  initially  5mm  thick  will 
compress  to  as  small  as  10pm.  This  is  smaller  than  the  typical  explosive 
particle  size  and  than  the  surface  irregularities  in  the  explosive  sample. 
On  the  scale  of  its  final  thickness,  then,  the  gap  may  not  be  regarded 
as  planar.  During  the  final  portion  of  gap  closure  air  will  flow  from 
regions  in  which  the  explosive  protrudes  above  its  mean  surface  elevation 
into  surface  defects.  This  is  illustrated  schematically  in  Figure  10. 
Geometrically,  the  problem  retains  its  one-dimensional  character  only 
on  a  local  scale  in  cylinders  whose  dimensions  are  small  compared  to 
dimensions  over  which  significant  variation  in  the  surface  elevation 
occurs.  As  the  piston  is  closed  air  crosses  these  cylinder  boundaries. 
Thus,  the  mass  of  the  air  residing  in  the  final  local  gap  differs  from 
that  in  the  original  cylinder.  Attempts  to  model  the  effects  of  defects 
of  various  sizes  showed  that  the  defect  size  alone  was  not  an  important 
determinant  of  peak  interface  temperature  but  that  only  local  mass 
addition  was  effective  in  increasing  the  temperature.  This  is  related  to 
the  defect  size  in  some  unknown  manner.  Thus,  the  mass  addition  which 
occurs  at  surface  defects  may  be  introduced  into  the  one-dimensional 
model  by  simply  specifying  the  final  to  initial  mass  ratio  m^/iiu  and 

assuming  a  linear  variation  during  gap  closure.  These  results  are 

PISTON 


Figure  10.  Schematic  of  Convergence  Effect 
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summarized  in  Figure  11  for  a  =  a  =10  (since  in  the  absence  of 

max 

boundary  effects  convergence  alone  cannot  cause  the  interface  temperature 
to  exceed  the  adiabatic  compression  limit).  This  also  increases  the  tem¬ 
perature  prediction  enough  to  explain  ignition. 


nr^/m. 


Figure  11.  Effect  of  Convergence  on  Interface  Temperature 

4.  Combined  Effects.  A  single  case  was  considered  to  demonstrate  the 
combined  effects  of  concentrated  energy  transport  enhancement,  rapid 
pressurization  and  convergent  air  flow.  Relatively  mild  values  of  the 
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parameters  describing  these  effects  were  chosen  as  follows: 


a 

max 


2 


s  =  1 


^  =  2xl03  GPa/s 
dt 


m. 

1 


The  initial  conditions  are  those  corresponding  to  case  B.  The  resulting 
maximum  interface  temperature  was  668K. 


5.  Diesel ing.  If  fine  particles  of  explosive  protrude  or  are  entrained 
into  the  air,  the  temperature  at  the  surface  of  the  particles  will  exceed 
the  planar  interface  temperature  when  the  particles  are  sufficiently  small. 
The  Comp-B  heated  layer  thickness  given  by  equation  4  is  a  measure  of  the 
particle  size  required  to  yield  this  effect.  This  thickness  is  plotted  as 
a  function  of  time  in  Figure  12.  Times  range  between  10  ps  and  10  ms.  The 
maximum  time  for  events  observed  in  the  activator  is  approximately  1  ms 
which  corresponds  to  a  heated  layer  thickness  of  about  30  pm.  This  may 
be  interpreted  to  mean  that  particles  whose  radius  is  less  than  30  pm 
exhibit  a  surface  temperature  higher  than  that  associated  with  a  semi- 
infinite  layer.  It  is  reasonable  to  believe  that  particles  of  this  size 
may  be  present  and  to  conclude  that  dieseling  as  a  sensitizing  mechanism 
may  not  be  ruled  out. 


IV.  CONCLUSION 

Two  analytical  models  to  predict  the  temperature  history  of  an 
explosive-air  interface  under  conditions  of  compressive  heating  have  been 
reported  herein.  The  adiabatic  compression  model  is  inadequate  to 
describe  the  activator  experiments  since  it  excludes  effects  of  press¬ 
urization  rate,  gap  thickness,  piston  properties  and  mass  addition  or 
leakage,  and  improperly  describes  the  effect  of  initial  pressure.  The 
finite  compression-rate  model,  on  the  other  hand,  yields  qualitative 
agreement  with  the  activator  experiments  with  respect  to  all  of  the 
above  stimulus  variations.  In  addition,  computations  for  the  alternate 
gases  -  argon,  carbon  dioxide  and  hydrogon  -  show  similar  effects  on 
sensitivity  using  both  models.  The  most  important  observation,  however, 
is  that  this  planar  one-dimensional  model  does  not  explain  ignition. 

This  may  indicate  that  multidimensional  processes  are  important  in  the 
experiments.  Several  such  processes  have  been  identified  and  introduced 
into  the  models.  The  most  important  of  these  is  turbulence  in  the 
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THICKNESS  OF  COMP  B  HEATED  LAYER  (  m ) 


r 


Figure  12.  Comp-B  Heated  Layer  Thickness 


air  which  serves  to  enhance  the  energy  transport  to  the  interface.  This 
effect  is  required  to  increase  the  theoretical  maximum  temperature  which 
is  otherwise  too  low.  The  finite  compress ion- rate  model  results  show 
that  enhanced  energy  transport  alone  still  does  not  explain  ignition 
unless  the  effect  is  concentrated  near  the  explosive-air  interface.  A 
combination  of  this  effect  with  either  rapid  pressurization  or  convergent 
airflow  is  necessary  when  no  such  boundary  effects  exist.  A  combination 
of  these  effects  will  produce  temperatures  high  enough  to  explain  ignition 
with  relatively  mild  deviation  of  each  influencing  factor  from  the  one¬ 
dimensional  planar  value.  Finally,  the  dieseling  mechanism  of  sensitization 
has  been  discussed.  The  particle  size  required  for  entrained  particles 
to  show  higher  surface  temperatures  in  the  time  frame  of  events  in  the 
activator  is  within  the  range  of  those  which  are  actually  present  in  the 
explosive  charge.  These  conclusions  are  subject  to  the  limitations  of 
the  model. 
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APPENDIX 


The  flow  of  the  gas  during  gap  closure  should  satisfy  both  the 
equations  of  continuity  and  conservation  of  momentum.  In  nondimensional 
form  these  are 


3p  _  +  XF  "U 

3t  F 


3fi._ 

3x 


F 


3u 

3x 


3u 

TT 


+  xF  -u 


3u  4 

*  *  Z? 


PF  3x 


However,  they  are  supplanted  in  the  present  model  by  the  assumptions  of 
linear  velocity  variation  and  uniform  pressure.  It  is  of  interest  to 
examine  the  conditions  under  which  these  assumptions  satisfy  continuity 
and  conservation  of  momentum. 

The  linear  velocity  variation  is  written 

u  =  xUp  =  U£  +  xF 

so  that 

+  xF  -u  ■  0  , 


and 
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since  =  0  for  the  gap.  With  these  expressions,  the  continuity 
equation  becomes 

i  ap  .  i 

p  St  F 

Since  the  right  hand  side  is  a  function  of  t  only,  density  variations 
with  x  must  be  negligible  and 

dp  -dF 

P  *  r~ 

or  the  density  varies  inversely  as  the  gap  thickness. 

Since  pressure  is  uniform,  =  0  and  the  momentum  equation  may  be 
• ....  oX 

written 

xF  =  0 


This  implies  that  the  piston  acceleration  must  be  negligible. 


NOMENCLATURE 


specific  heat  at  constant  pressure 

nondimensional  layer  thickness 
grid  point  index 

heat  produced  by  chemical  reaction 
reference  length 
mass  per  unit  area 
pressure 

heat  flux  per  unit  area 
gas  constant 
universal  gas  constant 
temperature 
time 

boundary  velocity 
particle  velocity 
distance 

heating  efficency 
layer  thickness 

ratio  of  specific  heats,  c^  /  cv 

thermal  conductivity 

density 

thermal  conductivity  multiplier 


viscosity 


SUBSCRIPTS 


air  value  of  constant  for  air 

f  value  after  piston  closure 

he  value  of  constant  for  explosive 

I  value  at  explosive  air  interface 

i  value  at  initial  state 

£  left  end  of  region 

o  reference  value 

p  value  for  piston 

r  right  end  of  region 

<*>  theoretical  maximum  value 


SUPERSCRIPTS 

k  previous  time  step 

£  present  time  step 

n  region  index 

—  dimensional  parameter 

.  derivative  with  respect  to  t 
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